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Dilepton production in In+In collisions at 158 A-GeV is studied within the microscopic parton- 
hadron-string dynamics (PHSD) transport approach that incorporates explicit partonic degrees-of- 
freedom, dynamical hadronization as well as the more familiar hadronic dynamics in the final reac- 
tion stages. A comparison to the data of the NA60 Collaboration shows that the measured dilepton 
yield is well described by including the coUisional broadening of vector mesons, while simultaneously 
accounting for the electromagnetic radiation of the strongly coupled quark-gluon plasma (sQGP) 
via off-shell quark-antiquark annihilation, quark annihilation with additional gluon Bremsstrahlung 
and the gluon-Compton scattering mechanisms. In particular, the spectra in the intermediate mass 
range (1 GeV < M < 2.5 GeV) are dominated by quark-antiquark annihilation in the nonpertur- 
bative QGP. Also, the observed softening of the transverse mass spectra at intermediate masses 
(1 GeV < M < 2.5 GeV) is approximately reproduced. Furthermore, for dileptons of low masses 
(M < 0.6 GeV), we find a sizable contribution from the quark annihilation with additional gluon 
bremsstrahlung, thus providing another possible window for probing the properties of the sQGP. 

PACS numbers: 25.75.-q, 13.60.Lc, 12.38.Mh, 14.40.Lbp, 14.65. Dw 



I. INTRODUCTION 

Electromagnetic probes, i.e. dileptons and photons, 
are powerful tools to explore the early hot, dense stage 
of heavy-ion collisions as they are essentially unaffected 
by final-state interactions. Through their invariant mass 
and momentum distributions, they carry to the detector 
information about the conditions and properties of the 
environment in which they are emitted, thus providing 
a glimpse deep into the bulk of the strongly interacting 
matter created in these collisions [l|, 13- In particular, 
dileptons have been suggested as probes of the quark- 
gluon plasma (QGP) that is expected to be produced 
during the early stage of heavy-ion collisions at Super- 
Proton- Synchrotron (SPS) energies Q-Q]- 

Recently, the NA60 Gollaboration Q has measured 
dileptons from In-|-In collisions at 158 A-GeV and found 
that the inverse slope parameter or effective temperature 
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of the transverse mass spectrum of dileptons in the in- 
termediate mass region is lower than that of dileptons at 
lower masses, which are dominantly of hadronic origin. 
This might be explained if the dilepton spectrum at in- 
variant masses above 1 GeV is essentially due to partonic 
channels in the QGP j9l-KLl|. In this case, the softening of 
the transverse mass spectrum with increasing invariant 
mass implies that the partonic channels occur dominantly 
before the collective radial flow has developed. 

Since dileptons are emitted over the entire history of 
the heavy-ion collision, - from the initial nucleon-nucleon 
collisions through the hot and dense phase and to the 
hadron decays after freeze-out, - microscopic covariant 
transport models are very useful for disentangling the 
various sources that contribute to the final dilepton spec- 
tra seen in experiments. The assumption that the dilep- 
ton spectrum at masses above 1 GeV might be domi- 
nated by radiations from the QGP was supported by 
studies within the Hadron-String-Dynamics (HSD) trans- 
port approach [l^l, which has shown [13| that the mea- 
sured dilepton yield at low masses (M < 1 GeV) can 
be well explained by dilepton production from hadronic 



interactions and decays, while there is a discrepancy be- 
tween the HSD results and the data in the mass region 
above 1 GeV. The excess seen for M > 1 GeV could 
not be accounted for by hadronic sources in HSD with 
or without medium effects and might be interpreted as 
a signal for the existence of partonic matter already in 
heavy-ion collisions at 158 A-GeV incident energy. In- 
deed, results from model studies by Dusling and Zahed 
[ij] as well as Renk and Ruppert [15| have indicated that 
this excess could be due to partonic channels, i.e. primar- 
ily to qq annihilation. On the other hand, this dilcpton 
excess has been attributed by Rapp and collaborators 
to multi-meson production channels (denoted shortly as 
47r-contribution [la|)- These different interpretations of 
the experimental data are still being extensively debated. 
In this respect, the physics of dilcpton transverse mo- 
mentum spectra can be especially relevant jl7{ . Due to 
the noncquilibrium nature of heavy-ion collisions, a clar- 
ification within a transport approach that incorporates 
dilcpton production from the (non-equillibrium) partonic 
phase, hadronic decays and the microscopic secondary 
hadronic interactions - including the "47r" channels - 
thus appears appropriate. 

Another open question to be answered by the micro- 
scopic transport calculations is the existence of other 
'windows' in the phase space for observing dileptons from 
the quark-gluon plasma over the hadronic sources. It 
has been originally suggested that the substantial ther- 
mal yield from the deconfined phase existed in the invari- 
ant mass region between the and J/^ peaks Q, while 
the spectrum at lower masses was dominated by meson 
decays. On the other hand, the calculations [13, U^ 
of the thermal dilcpton yield from qq annihilation in a 
blast wave model in comparison to that from a hadronic 
cocktail, the Drell-Yan mechanism and the correlated 
semileptonic decays of open charm have found a possi- 
ble second region of the phase space for the observation 
of this thermal source at masses w 0.3 — 0.6 GeV and low 
transverse momentum. 

The Parton-Hadron-String Dynamics [13, IHI (PHSD) 
transport approach, which incorporates the relevant off- 
shell dynamics of vector mesons and the explicit partonic 
phase in the early hot and dense reaction region as well as 
the dynamics of hadronization, allows for a microscopic 
study of various dilcpton production channels in noncqui- 
librium matter. The PHSD off-shell transport approach 
is particularly suitable for this investigation, since it in- 
corporates various scenarios for the modification of vec- 
tor mesons in a hot and dense medium, seen experimen- 
tally in the enhanced production of lepton pairs in the 
invariant mass range 0.3 < 1\I < 0.7 GcV/c^. In the 
present work, we calculate dilepton production from the 
partonic and hadronic sources within PHSD by including 
the multi-meson channels and the partonic channels be- 
sides the usual hadron decay channels. By consistently 
treating in the same microscopic transport framework 
both partonic and hadronic phases of the collision sys- 
tem, we are aiming to determine the relative importance 



of different dilcpton production mechanisms and to point 
out the regions in phase space where partonic channels 
are dominant. 

The paper is organized as follows. In Sec. II, we give 
a brief description of the PHSD approach. We then de- 
scribe in Sec. HI the partonic sources of dilepton pro- 
duction incorporated in PHSD and in Sec. IV dilepton 
production by (in-medium) hadrons and in multi-meson 
processes. In Sec. V, we compare the results of the cal- 
culations to the available experimental data. Finally, the 
conclusions are given in Sec. VI. 



II. PHSD TRANSPORT APPROACH 

To address dilepton production in a hot and dense 
medium as created in heavy-ion collisions, we employ 
an up-to-date relativistic trans por t model, i.e. the Par- 
ton Hadron String Dynamics [13, [HI (PHSD). PHSD 
consistently describes the full evolution of a relativistic 
heavy-ion collision from the initial hard scatterings and 
string formation through the dynamical deconfinement 
phase transition to the quark-gluon plasma as well as 
hadronization and to the subsequent interactions in the 
hadronic phase. 

In the hadronic sector, PHSD is equivalent to the 
Hadron-String-Dynamics (HSD) transport approach |12l . 
l22i . i23i] that has been used for the description of pA and 
AA collisions from SIS to RHIC energies and has lead to 
a fair reproduction of measured hadron abundances, ra- 
pidity distributions and transverse momentum spectra. 
In particular, the HSD incorporates off-shell dynamics 
for vector mesons [2J| and a set of vector-meson spec- 
tral functions [23 that covers possible scenarios for their 
in-medium modification. 

The transition from the partonic to hadronic degrees 
of freedom is described by covariant transition rates for 
the fusion of quark-antiquark pairs to mesonic resonances 
or three quarks (antiquarks) to baryonic states, i.e. dy- 
namical hadronization. Note that due to the off-shell 
nature of partons on one hand and the resulting hadrons 
on the other, the hadronization process obeys all conser- 
vation laws (i.e. 4-momentum conservation, flavor cur- 
rent conservation) in each event, the detailed balance re- 
lations, and the increase in the total entropy S. The 
transport theoretical description of quarks and giuons in 
the PHSD is based on a dynamical quasiparticle model 
(DQPM) for partons matched to reproduce lattice QCD 
(IQCD) results for a quark-gluon plasma in thermody- 
namic equilibrium. The DQPM provides the mean-fields 
for gluons/quarks and their effective 2-body interactions 
for the implementation to PHSD. 

We briefiy recall the basic assumptions of the DQPM 
model (for details about the DQPM model and the off- 
shell transport we refer to Ref. [2Q|). Following Ref. J27| . 
the dynamical quasiparticle mass (for giuons and quarks) 
is assumed to be given by the thermal mass in the 
asymptotic high-momentum regime, which is proper- 
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FIG. 1: (color online) The number of q + q (a) and q + g (b) collisions, in which a dilepton pair is produced, in a central In + In 
reaction at an incident energy of 158 A-GeV versus the masses of the (quasiparticle) quark and antiquark as resulting from 
PHSD. 



tional to the temperature T and a running coupling 
g(T/Tc) (squared), for which the following parametriza- 
tion is used 



9\T/T,) 



487r2 



(lliV, - 2Nf) ln[A2(r/r, - T,/re)2] 



(1) 



Here Nc = 'i stands for the number of colors while Nf 
denotes the number of flavors. The parameters control- 
ling the infrared enhancement of the coupling A and Tg 
have been fitted in Ref. [21| to recent IQCD results for 
the entropy density s{T). An almost perfect reproduc- 
tion of the energy density e{T) and the pressure P{T) 
from IQCD is achieved as well (cf. [2ll|). 

In line with Rcf . [23| , the parton spectral functions are 
no longer (5— functions in the invariant mass squared but 
taken as 
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separately for quarks and gluons {j = q^q,g). With the 
convention E'Hp) — p^+Mj—jj, the parameters Mj and 
7j are directly related to the real and imaginary parts of 



the retarded self-energy, e.g. 11^ 



Mf 



2i'~fjUj. 



The width for gluons and quarks (for vanishing chem- 
ical potential /i^) is adopted in the form 
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where c = 14.4 (from Ref. [281) is related to a magnetic 
cut-off. 

We stress that a non-vanishing width 7 is the main 
difference between the DQPM and conventional quasi- 
particle models [2^. Its influence is essentially seen in 



correlation functions, e.g., in the stationary limit of the 
correlation in the off-diagonal elements of the energy- 
momentum tensor T*^' which defines the shear viscosity 
77 of the medium [28| . Here a sizeable width is mandatory 
to obtain a small ratio of the shear viscosity to entropy 
density 77/s, which results in a roughly hydrodynamical 
evolution of the partonic system in PHSD [301 ■ The fi- 
nite width leads to two-particle correlations, which are 
taken into account b y m eans of the generalized^ off-shell 



transport equations j24| that go beyond the mean field 
or Boltzmann approximation [26l |3l| . 
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The off-shell effect can be seen for example in Fig. 
where the number of the q + q and q + g collisions 
which a dilepton pair can be formed - is shown as a func- 
tion of the participating parton masses. The plots have 
been generated by a simulation in PHSD for a central 
In + In reaction at an incident energy of 158 A-GeV. 
The maximum of the distribution indicates the average 
pole mass of the quark/gluon, while the width correlates 
with the average width of their spectral function. The 
values for the masses and widths are in agreement with 
those from the DQPM fit to the lattice data for the tem- 
peratures in the range « 1 — 2 Tc- 

For an illustration of the quark and gluon interactions 
in a heavy- ion collision as generated in PHSD, we show in 
Fig. m the number oi q + q (solid line) and q + g (dashed 
line) collisions that can create dilepton pairs per event 
in a central In + In reaction at an incident energy of 
158 A-GeV versus the invariant energy ^/s of the elemen- 
tary partonic collision. One can see that the tails of the 
collision distributions calculated in the PHSD transport 
are almost exponential, thus close to thermal. On the 
other hand, the collisions at very low ^/s are suppressed. 
This 'threshold effect' is due to the finite masses of the 
dynamical quarks, antiquarks and gluons. Additionally, 
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FIG. 2: (color online) Number of parton collisions per event in 
a central In + In reaction at an incident energy 158 A-GeV 
versus the invariant energy y/s of the elementary partonic 
collision as simulated in PHSD. The number oi q + q collisions 
is given by the solid (red) line while that oi q-\- g collisions is 
given by the dashed (blue) line. 



one notices that the threshold is not sharp because of 
the rather broad spectral functions (and therefore broad 
mass distributions) of the colliding partons. 



III. PARTONIC SOURCES OF DILEPTONS IN 
PHSD 

In the scope of the one- and two-particle interactions, 
dilepton radiation by the constituents of the strongly in- 
teracting QGP proceeds via the elementary processes il- 
lustrated in Figs. [3] and |31 the basic Born g -f g annihi- 
lation mechanism, gluon Compton scattering (<? + g — > 
7* -|- g and 9 + 5 — > 7* + g), and quark -f- anti-quark an- 
nihilation with gluon Bremsstrahlung in the final state 
{q + q ^ g + 7*), virtual quark decay (g — > g + g -I- 7*) 
and virtual gluon decay (5 — > g + g + 7*) • In the on-shell 
approximation, one uses perturbative QCD cross sections 
for the processes listed above. However, in the strongly 
interacting QGP the gluon and quark propagators differ 
significantly from the non-interacting propagators. Ac- 
cordingly, the cross sections for dilepton production in 
the partonic channels have been calculated in Ref. [33 
in the DQPM model that had been fitted to lattice QCD 
results in thermal equilibrium before [28| . 

The importance of finite mass corrections to the per- 
turbative cross sections has been stressed in Ref. [33 ). 
It was shown that the finite quark and gluon masses 
can modify the magnitude as well as the M— and 
Pt— dependence of the cross sections of the processes in 
Fig. 1 compared to the perturbative results for massless 
partons (cf. Figs. 3 and 4 of Ref. [S^l). The modifica- 
tions are large at lower M^ and at the edges of the phase 



space. It was shown that the most prominent effect of 
the quark masses on the dimuon production cross sec- 
tions in the Born mechanism (g + g — > 7*) was a sharp 
threshold value for the invariant mass of the dilepton pair 
Mmin = mi+m2. On the other hand, the finite masses of 
the quark and antiquark produce additional higher-twist 
corrections to the cross section, which decrease with in- 
creasing M^, so that the off-shell cross sections approach 
the leading twist - on-shell - result in the limit of high 
dilepton masses. In Fig. 4 of Ref. [32[, an analogous 
comparison for the 2 — > 2 process g -I- g — t- 7* + g was 
shown by plotting the off-shell (i.e. with finite masses 
for the quarks and gluons) cross section for the quark 
annihilation with gluon bremsstrahlung in the final state 
at various values of the quark and gluon off-shellnesscs 
(masses) and the corresponding on-shell result. As found 
in Ref. [32|, the maximum pair mass shifts to a lower 
value (in order to produce a massive gluon in the final 
state). For the rest of the M values, the effect of the 
quark and gluon masses is about 50%. For m^/g — > 0, 
the cross section approaches the leading twist pQCD re- 
sult. 

The question of the effect of a finite parton width - 
which parametrizes the effect of their interaction rate and 
correlation, including multiple scattering - on dilepton 
rates in heavy-ion collisions was addressed in Ref. [32| 
by convoluting the off-shell cross sections with phe- 
nomcnological spectral functions A{mq) and A{mg) for 
the quarks and gluons in the quark-gluon plasma and 
with parton distributions in a heavy-ion collision simi- 
lar to those of Fig. [5] in the present paper. The finite 
width of the quasiparticles was found to have a sizable 
effect on the dilepton production rates. In particular, the 
threshold of the Drell-Yan contribution was washed out. 
Also, the shape and magnitude of the 2 — > 2 processes 
(g -|- g — > g + 7* and q + g —>■ g + 7*) were modified. 
One further observed that the contribution of the gluon 
Compton process g 4- g — )■ g -I- 7* to the rates was small 
compared to that of g -f g annihilations. 

In the present work, we implement the cross sections 
obtained in Ref. [33 into the PHSD transport approach 
in the following way: Whenever the quark-antiquark, 
quark-gluon and antiquark-gluon collisions occur in the 
course of the Monte-Carlo simulation of the partonic 
phase in PHSD, a dilepton pair can be produced accord- 
ing to the off-shell cross sections [32[, which depend, in 
addition to the virtualities of the partons involved, on 
the energy density in the local cell, in which the collision 
takes place. The local energy density governs the widths 
of the quark and gluon spectral functions as well as the 
strong coupling (cf. Eqs. ^ and ([3|) that depend on tem- 
perature T which in turn is uniquely related to the energy 
density by the lattice QCD equation of state). Numeri- 
cally, one finds from a PHSD simulation of a heavy-ion 
collision at SPS energies that the running coupling as 
in the partonic phase is often of order 0(1) and thus the 
contribution of the higher-order Bremsstrahlung diagram 
is compatible in magnitude to the Born term. 
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FIG. 3: Diagrams contributing to dilepton production from the QGP: (a) Drell-Yan mechanism, (b) gluon-Compton scattering 
(GCS), (c) vertex correction, (d) gluon Bremsstrahlung (NLODY), where virtual photons (wavy lines) split into lepton pairs, 
spiral lines denote gluons, and arrows denote quarks. In each diagram the time runs from left to right. 
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FIG. 4: Diagrams contributing to dilepton production by virtual quasi-particles in addition to those presented in Fig. [S] Left: 
the decay of a virtual quark; Right; the decay of a virtual gluon. Virtual photons (wavy lines) split into lepton pairs, spiral 
lines denote gluons, and arrows denote quarks. 



IV. HADRONIC SOURCES OF DILEPTONS IN 
PHSD 



In the hadronic sector, PHSD is equivalent to the 
Hadron-String-Dynamics (HSD) transport approach [IJ, 
[22l . [23j . The implementation of the hadronic decays 
into dileptons (tt-, 77-, ?/-, lu-, A-, ai-Dalitz, p -^ l^l^, 
uj -^ 1+1-, (j) -^ 1+1-) in HSD (and PHSD) is described 
in detail in Refs. [IJ, [23. For the treatment of the lep- 
tonic decays of open charm mesons and charmonia we 
refer to Refs. [SJ, [S^- In the present paper, we extend 
the hadronic sources for dilepton production to include 
secondary multi-meson interactions by incorporating the 
channels ttw — >■ Z+Z~, ttci — )■ I'^l" , and pp — > /+/". 

The dilepton production by a (baryonic or mesonic) 
resonance R decay in HSD and PHSD can be schemati- 
cally presented in the following way: 



BB -^ RX 




(4) 


mB -^ RX 
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R^e+e-X, 




(6) 


R -^ mX, m - 


> e+e^X, 
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R -^ R'X, R' - 


-^ e+e-X, 


(8) 



i.e. in a first step a resonance R might be produced 
in baryon-baryon {BB) or meson-baryon (jnB) colli- 
sions dU, (O. Then this resonance can couple to dilep- 
tons directly ([5]) (e.g., Dalitz decay of the A resonance: 
A — >■ e^e~N) or decays to a meson m {+ baryon) 
or in (O produce dileptons via direct decays (p, w) or 



Dalitz decays (tt", 77,0;). The resonance R might also de- 
cay into another resonance R' ^ which later produces 
dileptons via Dalitz decay. Note, that in the combined 
model the final particles - which couple to dileptons - 
can be produced also via non-resonant mechanisms, i.e. 
'background' channels at low and intermediate energies 
or string decay at high energies. 



A. In-medium modification of vector mesons 

While the properties of hadrons arc rather well known 
in free space (embedded in the nonperturbative QCD 
vacuum) , the masses and lifetimes of hadrons in a bary- 
onic and/or mesonic environment are subject of current 
research that aims at achieving a better understanding 
of the strong interaction and the nature of confinement. 
For example, a broadening of the vector mesons can be 
understood as a shortening of the lifetime of the vec- 
tor mesons p, ui and in the medium. In this context 
the modification of hadron properties in nuclear mat- 
ter are of fundamental interest (cf. J35l - l39t ). since QCD 
sum rules [3^, 123, |40| as well as QCD inspired effec- 
tive Lagrangian models [35|, |3^, |4l|, 14^ predict significant 
changes, e.g., in the properties of the vector mesons (p, u 
and (j)) with the nuclear density pN and/or temperature 
T [lllil-lil. 

A modification of the properties of vector mesons in 
the nuclear medium was first seen experimentally in the 
enhanced production of lepton pairs above known sources 
in nucleus- nucleus collisions at SPS energies [4^, |43]. As 



proposed in Refs. [481, |49|, the observed enhancement 
in the invariant mass range 0.3 < M < 0.7 GeV/c^ 
might be due to a shift of the p-mcson mass following 
the Brown/Rho scaling [35[ or the Hatsuda and Lee sum 
rule prediction 36l. The microscopic transport studies in 
Refs. [IJ, [22I, [50, IMl for these systems have given support 
for this interpretation. On the other hand, more conven- 
tional approaches that describe a melting of the p-meson 
in the medium due to the strong hadronic coupling (along 
the lines of Refs. [4l|, |42|) have also been found to be 
compatible with the early CERES data [H, [IS [H, [H . 
This ambiguous situation has been clarified to some ex- 
tent in 2006 by the NA6G Collaboration since the invari- 
ant mass spectra for /i^/-t^ pairs from In-|-In collisions 
at 158 A-GeV favored the 'melting p' scenario [g. Also, 
the later data from the CERES Collaboration (with en- 
hanced mass resolution) [5J] showed a preference for the 
'melting p' picture. 

The various models, which predict a change of the 
hadronic spectral functions in the (hot and dense) nu- 
clear medium, may be classified into two different cate- 
gories: i) a broadening of the spectral function or ii) a 
mass shift of the vector mesons with density and/or tem- 
perature. In view of many-body dynamics, both modifi- 
cations should be studied simultaneously as well [l|, [531 ■ 
Thus we explore in the present study three possible sce- 
narios with respect to the low-mass dilepton spectrum 
which essentially addresses all possible properties of the 
p-meson: (1) a broadening of the p spectral function, (2) 
a mass shift, and (3) a broadening plus a mass shift. 

We incorporate the effect of coUisional broadening of 
the vector-meson spectral functions (as in Refs. [56l[57|). 
by using for the vector meson width 



r*y{M, \p\,pn) = rv{M)+T,oii{M, \p\,pn)- 



(9) 



Here ry(M) is the total width of the vector mesons 
{V = p, to) in the vacuum. The coUisional width in ([9]) is 
approximated as 



^coIi{M,\p\,pn) 



-f PN <v a*^l; >? 



acoii 



Pn 



(10) 



Here v = \p\/ E] p, E are the velocity, 3-monientum and 
energy of the vector meson in the rest frame of the nu- 
cleon current and 7^ ~ 1/(1 — w^). Furthermore, pN 
is the nuclear density and (Ty% the meson-nucleon total 
cross section. 

In order to simplify the actual calculations for dilep- 
ton production, the coefficient Ucoii has been extracted in 
the PHSD transport calculations from the vector-meson 
collision rate in In + In reactions at 158 A- GeV as a 
function of the density p^- In case of the p meson the 
collision rate is dominated by the absorption channels 
pN — > -kN or pN — > Att — > mrN . Also the reactions 
p + TT <-> oi are incorporated. The numerical results for 
^coIi{pn) then have been divided by pn / Po to fix the co- 
efficient acoii h^ (fTO]). We obtain acoii « 150 MeV for the 
p and acoii ~ 70 McV for uj mesons which are consistent 



with Ref. [58J. In this way the average effects of col- 
lisional broadening arc incorporated in accordance with 
the transport calculations and allow for an explicit rep- 
resentation of the vector-meson spectral functions versus 
the nuclear density as demonstrated in Ref. [ZSj . 

In order to explore the observable consequences of 
vector meson mass shifts at finite nuclear density, the 
in-medium vector meson pole masses are modeled (op- 
tionally) according to the Hatsuda and Lee [3g or 
Brown/Rho scaling J35| as 



m;{p 



Mo 



N) 



(1 + apN/po) ' 



(11) 



where pN is the nuclear density at the resonance decay 
position r; po = 0.16 fm~ is the normal nuclear density 
and a ~ 0.16 for the p and a ~ 0.12 for the a; meson 
[SSj . The parametrization pT|) may be employed also 
at much higher collision energies and one does not have 
to introduce a cut-off density in order to avoid negative 
pole masses. Note that (fTT|) is uniquely fixed by the 
'customary' expression Mq{pn) w il/o(l — apN/po) in 
the low density regime. 

The spectral function of the vector meson V for the 
mass M at baryon density pN is taken in the Breit- 
Wigner form: 



Ay(M,pjv) = Ci 



APT*yiM,pN) 



-7T(M^-MS\pN)r + {Mr*yiM,pNW 

(12) 
The factor Ci is fixed by the normalization condition for 
arbitrary piy: 



Ml, 



Av{M,pn) dM = 1, 



(13) 



Mr, 



where Mum = 2 GeV is chosen as an upper limit for 
the numerical integration. The lower limit for the vac- 
uum spectral function corresponds to the two-pion de- 
cay, Mmin = 2TO,r, whercas for the in- medium coUisional 
broadening case Mmin = ^m,, — > with nie denoting the 
electron mass. Mg is the pole mass of the vector meson 
spectral function which is Mq (pn = 0) = Mg in vacuum, 
however, may be shifted in the medium for the dropping 
mass scenario according to Eq. (fTTj) . The resulting spec- 
tral functions for the p and a; meson are displayed in 



Fig. 2 of Ref. [25 1 . 

With increasing nuclear density pN elastic and inelas- 
tic interactions of the vector mesons shift strength to low 
invariant masses. In the 'coUisional broadening' scenario 
we find a dominant enhancement of strength below the 
pole mass for the p meson while the lo meson spectral 
function is drastically enhanced in the low- and high- 
mass region with density (on expense of the pole-mass 
regime). In the 'dropping mass -I- coUisional broaden- 
ing' scenario both vector mesons dominantly show a shift 
of strength to low invariant masses with increasing p^. 
Qualitatively similar pictures are obtained for the (j) me- 
son but quantitatively smaller effects are seen due to the 



lower effect of mass shifts and a substantially reduced 
(j)N cross section which is a consequence of the ss sub- 
structure of the (j) meson. 

Note that, just as the HSD, the PHSD incorporates 
the off-shell propagation for vector mesons - according 
to Ref. 12J]. In the off-shell transport, the hadrou spec- 
tral functions change dynamically during the propaga- 
tion through the medium and evolve towards the on-shell 
spectral functions in the vacuum. The PHSD off-shell 
transport approach is particularly suitable for investigat- 
ing the different scenarios for the modification of vector 
mesons in a hot and dense medium. As demonstrated 
in Ref. [25|, the off-shell dynamics is important for res- 
onances with a rather long lifetime in the vacuum but 
strongly decreasing lifetime in the nuclear medium (es- 
pecially w and (j) mesons) and also proves vital for the 
correct description of dilepton decays of p mesons with 
masses close to the two pion decay threshold. For a de- 
tailed description of the off-shell dynamics we refer the 
reader to Refs. [H, H El HI ■ 



where 
Pcm = \/{s' 



(rni + m2)2)(s - (toi — TO2)^)/2-v/s (15) 



is the center-of-mass momentum_of the colliding hadrons 
with the masses mi, m2, and |A^p can be written as 



\M\'^^ -T 



Ana 
g2 



^fJ.L>J^ I 



(16) 



with q = Pi +P2 = P3 +P4 ^nd the fine structure constant 
a. In (|16p . Lf^^ is the leptonic tensor given by 



^t.^ 



P3P4 



PiP3 



g^^P3-p, + mf), (17) 



while H'^'^ is a hadronic tensor for the reaction. 
The hadronic tensor H'^'^ for the reaction 7r+ 
e+ + e~ is given by 



ff^'' = (p^-p^)(p^-K; 



(18) 



which leads to the well-known result for the tttt annihila- 
tion cross section 



B. Multi-meson channels of dilepton production 

The dilepton excess yield in In + In collisions at 
160 A-GeV incident energy for M > 1 GeV/c^ was found 
to be dominated by partonic sources within the dynam- 
ical studies of Renk and Ruppert [l3 as well as Dusling 
and Zahed [1J|- On the other hand, the model of van 
Hees and Rapp [lal suggests a dominance of hadronic 
sources dubbed '47r channels'. In order to clarify this 
question, we have incorporated in the PHSD the '47r 
channels' for dilepton production on a microscopic level 
rather than assuming thermal dilepton production and 
incorporating a parametrization for the inverse reaction 
fj,'^ + uT — ^ Air's by employing detailed balance as in 
Refs. U^^. 

By studying the electromagnetic emissivity (in the 
dilepton channel) of the hot hadron gas, it was shown 
in Refs. [60J, |6l| that the dominating hadronic reactions 
contributing to the dilepton yield at the invariant masses 
above the peak are the two-body reactions, i.e. tt -(- p, 
TT + uj, p + p, TT + ai. This conclusion was supported by 
the subsequent study in a hadronic relativistic transport 
model |62l |. Therefore, we implement the above listed 
two-meson dilepton production channels in the PHSD ap- 
proach. In addition, some higher vector mesons (p' etc.) 
are tacitly included by using phenomenological form fac- 
tors adjusted to data. 

We determine the cross sections for the mesonic in- 
teractions with dileptons in the final state using an 
effective Lagrangian approach, following the works of 
Refs. [60, |62[- The dilepton production cross section is 
given by the product of a form factor and the square of 
a scattering amplitude 
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where M is the mass of the lepton pair, and mi is the 
mass of the lepton. The electromagnetic form factor 
|F^(Af)p plays an important role in this process, pro- 
viding empirical support for the vector meson dominance: 
the pion electromagnetic form factor is dominated by the 
p(770) meson. In Ref. [6j], Gale and Kapusta proposed 
the form, 



\FAM)\' = 



(20) 



(A/2_m;2)2 + (rn,r,)2' 

where m^ = 0.775 GeV, m'^ ^ 0.761 GeV, and P^ = 
0.118 GeV. 

According to Ref. |6J| , the cross section for irp annihi- 
lation is given by 



o-(7r+/9 -^ II) = 
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(21) 

Note that the cross section (PT|) is evaluated in the 
narrow-width approximation for illustration purposes 
only. This simplification is not used in the actual 
transport calculation. The electromagnetic form factor 
|F7rp(M)p can then be determined by analyzing the ex- 
perimental data for e~^e~ — >■ tt^tt~tt^. In Ref. |64{ . three 
isoscalar vector mesons, (f>{lQ2Q), w(1420), and a;(1670) 
were found to be important in order to fit the experimen- 
tal data, namely. 



F^piM) 



f gvTTf. 
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(22) 



Here the summation runs over the three vector mesons 
listed above. While the coupling constants g^ and g^Trp 
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FIG. 5: (color online) a) Cross sections for the reactions 
e"^e~ — >■ TT + p (a) and e^e^ -^ tv + uj (b) in our model versus 
the experimental data. Panel (c) gives the measured cross 
section of the e^e~ — >■ tt^tv'tt tv reaction versus the sum of 
the model cross sections for e'^e~ -^ -k + u, e^e~ — > tt + ai 
and e^e~ -^ p + p. 



can be determined from the measured widths, the cou- 
pling constants for other two mesons and the relative 
phases were determined by a fit to the experimental data 
ofRefs. [65','6(|. These coupling constants were extracted 
from the latest data of the DM2 collaboration [631 ^^^ 
the ND collaboration [g^. The parameters are listed in 
Ref. [6J] . The comparison of the fit to the experimental 



data is shown in Fig. 5(a 



The cross section for lepton pai r production in pion- 
omcga annihilation is given by |62| 



a{7r°Lu -^ II) = 
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(23) 

The form factor can be parametrized in terms of three 
isovector p-like vector mesons, p(770), p(1450), and 
p(1700). 
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Here the summation runs over the three p-like resonances 
listed above. The parameters used are tti^i = 0.77 GeV, 
mr2 = 1.45 GeV, m^g = 1.7 GeV, T^i = 0.118 GeV, 
Tr2 = 0.25 GeV, T^s = 0.22 GeV, Ai = 0.85, A^2 = 
—0.077, Ar3 — 0.034, where Ay ~ (gy7rw/5v)exp{i(/)y}. 
The comparison with the experimental data of the 
ND [Hi and ARGUS collaborations [gi is shown in 



Fig. 5(b) 



Additionally, we consider the reactions ttoi — J^ II and 
pp -^ II, which are effectively four-pion processes. Using 
the Lagrangian for the irai interaction 
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one obtains for the cross section of the ttai — > II process. 
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(26) 



where the value of the coupling constant g = {gp/ fp) is 
adjusted so that the experimentally measured radiative 
decay widths are reproduced. 

We obtain the hadronic tensor H^^^ for the reaction 
p+p~ — 7> e^e~ by generalizing the formula of Ref. [601 ^^ 



explicitly take into account the broad spectral functions 
of the colliding p-mesons: 



+h';^^p,0P2ah';'''p2,Pis/{m%m%), (27) 



with 

/^r^ = (p^-pOff"'' + (9"-P2)/'^ + (pf-/).9'^"- (28) 

In this case, the hadronic tensor depends on (generally 
different) masses of the colliding particles nipi and mp2 . 
In the actual transport calculations, nipi are distributed 
according to the dynamical spectral functions. Using 
(|T4l) - (|16l) and ((27)) we obtain the following cross section 
as a function of M, nipi and mp2 



na'\Fpp\' 



aipp — >■ I'^l ) = „ „ ^ 

^' ^ ' l20ml,m%M^pcm 

X {9m^i + 18TO^i(3TOp2 - 2A/2) 

+im-'p2 - M^ [819to^2 + GS^ml^M"^ - IIM"*] 
-2ml^ [363m^2 + 32to^2M^ + 327ml^M^ - 2M^] 
+m^i [-156m^2 + 266ml2M^ + MM^] } , (29) 

which reduces in the narrow width approximation to 



a{pp -> I'^l ) 



na'lFppl^ 



60m4M3,/M2-4m2 



X (840m^ + 1076TO^Af2 - ebSm^M^ + 11 Af^). (30) 

The form factors |i^,rai P and l-FppP can be determined 
by analyzing the e"*'e^ — > 7r+7r~7r+7r^ and e'^e~ —^ 
7r"'"7r~7r'^7r° data. We determine |F7rai(Af)p from exper- 
imental data for e+e~ ^ 7r+7r^7r+7r~ from the 772 col- 
laboration Izi , the MSN collaboration fzl , and the ND 



collaboration [68|. Further constraints on |F^ai P and the 
determination of |i^ppP were provided by the experimen- 
tal data for e~^e~ — >■ 7r+7r~7r'^7r'', which can come from 
TTw, TTtti and pp intermediate states. Our form factors 
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(31) 



with rriri = 0.77 GeV, mr2 = 1-45 GcV, m^a = 1.7 GeV, 
Tri = 0.118 GcV, Tr2 == 0.25 GcV, Trz = 0.235 GeV, 
Ari = 0.05, Ar2 = 0.58, Ars. = 0.027, and 



Fpp{M) = Y. 



^i4>Vrn? 
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\ 9v J (m^ - M2) - imvTv 



V 



(32) 



with rriri = 0.77 GeV, 771^2 = 1.45 GeV, m^s = 1.7 GeV, 
r^i = 0.118 GeV, r^2 = 0.237 GeV, T^g = 0.235 GeV, 
Ari = 0.05, Ar2 = 0.05, 4r3= 0.02. The comparison to 



the data is shown in Fig. 5(c) 



Let us summarize that in order to fix the form fac- 
tors in the cross sections for dilepton production by the 
interaction of tt + p, n + co, p + p and ttai, we use the 
measurements in the detailed-balance related channels: 
e+e" — > TT -I- p, e~^e~ — >■ tt -I- a;, e+e" -^ p + p, and 
e+e^ ^^ TT+ai . Note that we fitted the form factors while 
taking into account the widths of the p and ai mesons in 
the final state by convoluting the cross sections with the 
(vacuum) spectral functions of these mesons in line with 
Ref. [721 (using the parametrizations of the spectral func- 
tions as implemented in HSD and described in [73|). In 
Fig. [5] we present the resulting cross sections, which are 
related by detailed balance to the ones we implemented 
into PHSD. 



V. DILEPTON RATES IN THERMAL 
EQUILIBRIUM 

Before proceeding to the results of the transport cal- 
culations and the comparison to data, we dedicate this 
section to a study of the dilepton spectrum qualitatively 
in a thermal model. In Fig. [51 the dilepton production 
rates in thermal equilibrium are presented. We assume 
here that the system evolves through a thermalized sys- 
tem of quark in the hot initial stage of the heavy-ion 
collision and through the state of a high-density hadron 
gas in the later phase of the collision. 

The main elementary process of dilepton production 
in a hadron gas is the pion annihilation into dileptons, 
mediated through vector meson dominance by the rho 
meson (tt + tt ^ p ^ 7* — > Z"*" -|- /~) and controlled by 
the pole at the rho mass of the pion electromagnetic form 
factor. For the pion annihilation, we use the standard 
cross section as, e.g., in Ref. 160| and the Breit-Wigner 
form factor with the pole mass and width of the p meson. 

In the partonic sector, the main sources of the dilep- 
tons are the reactions of quark-antiquark annihilation 
with the production of the virtual photon. Consider- 
ing the temperatures and baryon densities relevant for 
the SPS energies, we expect in PHSD the contribution 
of the processes involving gluons to be small compared 
to the leading q + q mechanism of dilepton production 
(note, however, that at higher energies, such as those of 
RHIC and LHC, gluons can play an important role in 
the dilepton production [7J). For the calculation of the 
QGP yield in the qualitative analysis of this section, the 
most simple perturbativc QCD cross sections are used 
for the processes q + q ^ l^l'' and q + q ^ g + l^l^, 
assuming as = 0.8. 

Thus wc plot the dilepton yields from the reactions 
TT + TT and q + q, where the pious and quarks have in 
general different temperatures T^ and Tqqp. The space- 
time volumes of the two phases are assumed to be approx- 
imately equal. In Fig. 6(a)[ the gas of pious is assumed 
to have the temperature T^ = 150 MeV, while the gas of 
quarks the temperature Tqgp = 250 MeV. In Fig. |6(b)[ 
we have Tqgp = T^ = 190 MeV. 
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FIG. 6: (color online) Rates of dileptons created in q + q and tt + ir annihilations within a thermalized gas of quarks with 
temperature Tqgp and pious with temperature T^. (a) Tqgp = 250 MeV, T^ = 150 MeV; (b) Tqgp = T,r = 190 MeV. 



It has been originally suggested that a 'window' for ob- 
serving dileptons from the plasma exists in the invariant 
mass region between the 4> and J/^P peaks Q. This is 
supported by the results shown in both Figs. 6(a) and 
|6(b)[ However, we see in Fig. 6(a) another region, i.e. 
M < 0.5 GeV, in which the q + q annihilation is com- 
patible or even larger than the radiation from the tt + tt 
annihilation; the contribution of the two-to-two process 
q + q ^- g + l~^l~ is especially important. The dominance 
of the thermal yield from quark interactions at masses 
below sa 0.5 GeV is in agreement with the conclusions 
of [13, 13 • The transport model results of the next sec- 
tion will clarify which of the equilibrium scenarios pre- 
sented in Fig. |6]- (a) or (b) - gives a closer resemblance 
to the channel decomposition of the dilepton production 
within a microscopic simulation. 

Of course, the observation of the QGP channels at low 
mass is possible only after the dilepton yield from the 
TT-, 77- and cj-Dalitz decays is removed. Another word of 
caution is in place here, because in Figs. [HI the vacuum 
properties of the p-meson have been used in plotting the 
TT -|- TT contribution, whereas the p-mcson properties are 
expected to be modified in medium. The modification of 
the p will change the size of the new, low mass window 
of the QGP observation (cf. next section). 



VI. RESULTS AND COMPARISON TO DATA 

Let us first note that the bulk properties of heavy-ion 
reactions at the top SPS energy, such as the number of 
charged particles, as well as their rapidity, pr and trans- 
verse energy distributions , were rather well described by 
PHSD; we refer to Ref. [20] for an extended and detailed 
comparison to the data. As the IQCD equation of state 



employed here has a crossover transition, the PHSD cal- 
culations show a rather long QGP phase in central In+In 
coUisions at 158 A-GeV (cf. Fig. 10 of Ref. [13) with the 
partonic degrees of freedom dominating for about 3 fm/c. 
Also, the elementary pp channel is well under control in 
PHSD as has been demonstrated in Ref. 



13|. 



Previously, by employing the HSD approach to the low 
mass dilepton production in relativistic heavy-ion colli- 
sions, it was shown in Ref. [ijl that the NA60 Collabo- 
ration data for the invariant mass spectra of ^"^p~ pairs 
from In-|-In collisions at 158 A- GeV favored the 'melting 
p' scenario @ . Also the data from the CERES Collabora- 
tion [5JI showed a preference for the 'melting p' picture. 
For other vector mesons {u>,(j)), the effects are relatively 
small, since, due to their much longer lifetimes, uj and (j) 
decay predominantly outside the medium after regaining 
the vacuum properties. 

As we see in Fig. [71 the current calculation in the PHSD 
approach confirms the earlier finding in the hadronic 
model that the NA60 data favor the scenario of the in- 
medium broadening of vector mesons. A comparison of 
the transport calculations to the data of the NA60 Col- 
laborations points towards a 'melting' of the p'-meson at 
high densities, i.e. a broadening of the vector meson's 
spectral function in line with the findings by Rapp and 
collaborators [7g . No pronounced mass shift of the vec- 
tor mesons is visible in the data. Thus the experimental 
results suggest that the approach to the chiral transition 
proceeds through broadening, and eventually melting, of 
the resonances rather than by dropping masses. On the 
other hand, a closer inspection of Fig. [3 shows that the 
conventional hadronic sources do not match the mea- 
sured yield at invariant masses above about 1 GeV/c^, 
while the yield at masses close to 1 GeV is reproduced 
by taking into account the dilepton production channels 
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FIG. 7: (color online) The HSD results for the mass differential dilepton spectra from In + In collisions at 158 A-GeV in 
comparison to the excess mass spectrum from NA60 [31 . The actual NA60 acceptance filter and mass resolution have been 
incorporated [7al • The solid lines show the HSD results for a scenario including the collisional broadening of the p-meson whereas 
the dashed lines correspond to calculations with 'free' p spectral functions for reference. The dash-dotted lines represent the 
HSD calculations for the 'dropping mass + collisional broadening' model. The (blue) bands represent the PHSD results 
incorporating direct dilepton radiation from the QGP in addition to a broadened p-meson. 



in the QGP. 

The NA60 collaboration has recently published accep- 
tance corrected data with subtracted charm contribu- 
tion [731 • In Fig. [5] we present PHSD results for the dilep- 
ton spectrum excess over the known hadronic sources as 
produced in In + In reactions at 158 A- GeV compared 
to the acceptance corrected data. We find here that the 
spectrum at invariant masses in the vicinity of the p peak 
is well reproduced by the p meson yield, if a broadening 
of the meson spectral function in the medium is assumed, 
while the partonic sources account for the yield at high 
masses. 

One concludes from Fig.|S]that the measured spectrum 
for M > 1 GeV is dominated by the partonic sources. In- 
deed, the domination of the radiation from the QGP over 
the hadronic sources in PHSD is related to a rather long 
- of the order or 3 fm/c - evolution in the partonic phase 
(in co-existence with the space-time separated hadronic 



phase) on one hand, cf. Fig. 10 of Ref. [201, and the 
rather high initial energy densities created in the colli- 
sion on the other hand, cf. Fig. 6 of Ref. |3J] . 

In addition, we find from Fig. [5] that in PHSD the 
partonic sources have a considerable (about 30%) contri- 
bution to the dilepton yield at M < 0.6 GeV. The yield 
from the two-to-two process q + q — >■ g + l~^l~ is especially 
important. This conclusion from the microscopic calcu- 
lation is in qualitative agreement with the conclusion of 
an early (more schematic) investigation in Ref. [T^j • 

Recalling the illustrative study of dilepton rates in 
thermal equilibrium in section |Vl we observe that the 
non-equilibrium microscopic simulation within the PHSD 
transport approach qualitatively implies a situation in 
which the initial partonic phase has temperatures of the 
order of Tqqp ~ 250 MeV and the hadron gas in the sub- 
sequent evolution a temperature Tjr ~150 MeV (assum- 
ing thermalization and that their evolution is approxi- 
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FIG. 8: (color online) Acceptance corrected mass spectra of 
excess dimuons from In + In at 158 AGeV integrated over 
Pt in 0.2 <Pt < 2.4 GeV from PHSD compared to the data 
of NA60 [731 . The dash-dotted line shows the dilepton yield 
from the in-medium p with a broadened spectral function, 
the dashed line presents the yield from the q + q annihila- 
tion, the dash-dot-dot line gives the contribution of the gluon 
Bremsstrahlung process (qq — >■ gl^l~), while the solid line 
is the sum of all contributions. For the description of the 
other lines, which correspond to the non-dominant channels, 
we refer to the figure legend. 



mately as long). A model scenario, in which the tempera- 
tures of the partonic and hadronic phases are equal for an 
extended period of space-time (T^r = Tqqp = 190 MeV) 
is not supported by the microscopic simulations. 

In order to elucidate the relative importance of the 
different hadronic sources of the excess dileptons in the 
heavy- ion collisions at top SPS energies, we show in Fig. [5] 
the channel decomposition of the main hadronic contri- 
butions to the dilepton rates in central In + In collisions 
at 158 A-GeV integrated over rapidity and pr- In partic- 
ular, the dilepton yield from the decays of the p-mesons 
(solid line) is dominantly composed of two channels: the 
direct p-mesons produced in mcsonic and baryonic strings 
(dashed line) and the 'thermal' p-mesons produced in 
TT -|- TT annihilations (dash-dotted line) . For comparison, 
the contributions of the '47r' processes are shown by the 
lines with symbols: the 7r+ai — > l^l^ process is displayed 
by the line with stars, tt + w — >■ Z+Z~ by the line with 
squares and p+p — )■ l'^l~ by the line with circles. Wc find 
that the dilepton yield from the decays of the 'thermal' 
p-mesons falls exponentially at high masses. The con- 
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FIG. 9: (color online) Dilepton radiation from p-mesons of 
different origins in PHSD from central In 4- In collisions at 
158 A-GeV compared to the contributions from the '47r' pro- 
cesses (ai -\--n, ■K-\-u}, p-\- p). The direct p-mesons produced in 
mesonic and baryonic strings are given by the dashed line and 
the 'thermal' p-mesons produced in 7r-|-7r annihilations by the 
dash-dotted line. The contributions of the '47r' processes are 
shown by the lines with symbols: the vr + ai — > l'^l~ process 
is displayed by the line with stars, tt + to — > I'^l^ by the line 
with squares and p + p ^>- l^l^ by the line with circles. 



tributions from the '47r' processes start dominating over 
the yield from the 'thermal' p decays at M « 1 GeV. We 
further confirm in PHSD that at i\/ > 1 GeV the con- 
tribution of the TT + ai process is the highest among the 
secondary mesonic sources of the dileptons, as was first 
noted by the authors of Ref. [60|. On the other hand, 
in contrast to the 'thermal' p's, the direct p-mesons pro- 
duced in the string decays (following the initial hard col- 
lisions) exhibit a power-law tail at masses above 1 GeV 
and, consequently, dominate the overall dilepton spec- 
trum of hadronic origin for M > 1 GeV. 

Next we investigate the centrality dependance of the 
dilepton production in heavy-ion collisions as SPS en- 
ergies. In Fig. [TO] we present the mass spectra of ex- 
cess dimuons from In + In at 158 AGeV for 0.2 < 
PT < 2.4 GeV and 3 < 77 < 4.2 from PHSD for differ- 
ent centrality bins. The dash-dotted, dashed, and solid 
lines show, respectively, the dilepton yield from the in- 
medium p with a broadened spectral function, the dilep- 
ton yield from the q + q annihilation and the sum of 
them. We have chosen the following centrality classes: 
central collisions (impact parameter 0.5 fm< b < 3.5 fm), 
semi-central (3.5 fm< b < 5.5 fm), and semi-peripheral 
(5.5 < b < 8.5 fm). The predictions in Fig. [TO] can be 
verified/falsified in the future by a direct comparison to 
the data as the latter become available. 

The yields in Fig. [TO] arc normalized to the number of 
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FIG. 11: (color online) Number of charged particles as a 
In+In, 158 A GeV, semi-peripheral (5.5<b<8.5) | function of the impact parameter from HSD (solid Unes) 

and PHSD (dash-dot lines) integrated over rapidity (upper 
lines) and within the pseudo-rapidity acceptance window of 
the NA60 experiment (lower lines). The vertical dashed lines 
indicate the different centrality classes. 
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FIG. 10: (color online) Mass spectra of excess dimuons from 
In + In at 158 AGeV for 0.2 < pr < 2.4 GeV and 3 < 77 < 4.2 
from PHSD for different centrality bins. The dash-dotted, 
dashed, and solid lines show, respectively, the dilepton yield 
from the in-medium p with a broadened spectral function, 
the dilepton yield from the q + q annihilation and the sum of 
them. 





< dNch/dr] > 


Centrality 


HSD 


PHSD 


b < 8.5 fm 


83.44 


79.00 


0.5 < 6 < 3.5 fm 


166.6 


157.1 


3.5 < 5 < 5.5 fm 


119.5 


112.6 


5.5 < 6 < 8.5 fm 


58.13 


55.54 



charged particles Nch- By studying the dependence of 
Nch on the centrality in In + In collisions at 158 AGeV 
in Fig. [TT] we find that PHSD and HSD give very simi- 
lar results (with only 5% quantitative difference). This 
finding is in line with the conclusions of the extended 
study in Rcf [20| that the multiplicities, rapidity- and 
transverse-momentum-distributions of the non-strange 
particles produced in heavy-ion collisions at 158 A-GeV 
are only weekly sensitive to the presence of a partonic 
phase with a (cross over) phase transition. The average 
numbers of charged particles per unit of pseudo-rapidity 
in PHSD and HSD for the chosen centrality classes are 
shown in the following table 1: 



Table 1 : The average numbers of charged particles per unit of 
pseudo-rapidity in PHSD and HSD for the different centrality 
classes. 

The NA60 Collaboration has accessed the information 
on the transverse momentum dependence of dilepton pro- 
duction by measuring the dilepton yield in different bins 
of pT- In Fig. [12] we show the mass spectra of excess 
dimuons from In + In at 158 AGeV for different trans- 
verse momentum bins from PHSD compared to the data 
of the NA60 Collaboration [1, ^. The dash-dotted, 
dashed/dash-dot-dot, and solid lines show, respectively, 
the dilepton yield from the in-medium p with a broad- 
ened spectral function, the dilepton yield from the q + q 
annihilation and the sum of them. One observes a gen- 
erally good agreement with the data. 

In Fig. I13[ transverse mass spectra of dilcptons for 
In+In at 158 AGeV in PHSD are compared to the data 
of the NA60 Collaboration for the four mass bins. The 
comparison of the mass dependance of the slope parame- 
ter evolution in PHSD and the data is shown explicitly in 
Fig. [M] Including partonic dilepton sources allows us to 
reproduce in PHSD the ?7i7--spectra (cf. Fig. ^1) as well 
as the finding of the NA60 Collaboration [^ [73] that the 
effective temperature of the dilcptons (slope parameters) 
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FIG. 12: (color online) Acceptance corrected mass spectra of excess dimuons from In + In collisions at 158 A-GeV for different 
transverse momentum bins from PHSD compared to the data of the NA60 Collaboration [3,[73- The dash-dotted, dashed/dash- 
dot-dot, and solid lines show, respectively, the dilepton yield from the in-medium p with a broadened spectral function, the 
dilepton yield from q + q and the sum of them. 



in the intermediate mass range is lower than that of the 
dileptons in the mass bin 0.6 < M < 1 GeV, which is 
dominated by hadronic sources (cf. Fig.[T4|). The soften- 
ing of the transverse mass spectrum with growing invari- 
ant mass implies that the partonic channels occur dom- 
inantly before the collective radial flow has developed. 
Also, the fact that the slope in the lowest mass bin and 
the highest one are approximately equal - both in the 
data and in PHSD - can be traced back to the two win- 
dows of the mass spectrum that in our picture arc influ- 
enced by the radiation from the sQGP: M = — 0.6 GcV 
and M > 1 GcV (cf. the discussions of Fig.[S]and Fig. [5]). 
A detailed look at the PHSD results shows that in total 
we still slightly underestimate the slope parameter Tg// 
in the p-mass region which might be due to missing par- 
tonic initial state effects or an underestimation of flow in 
the initial phase of the reaction. 



VII. SUMMARY 



In+In, 158GeV,LMR 



PHSD NA60 

■ 0.2<M<0.4, 24* 

• 0.4<M<0.6, 6x : 

• 0.6<M<0.9 

♦ 1.0<M<1.4 




0.4 0.8 1.2 

m -M [GeV] 



1.6 



FIG. 13: (color online) Transverse mass spectra of dileptons 
for In-|-In at 158 A-GeV in PHSD compared to the data of 
the NA60 Collaboration ^^. 



To address dilepton production in a hot and dense 
medium as created in heavy-ion collisions, we have em- 
ployed an up-to-date rclativistic transport model, the 
Parton-Hadron-String Dynamics ^^ (PHSD). PHSD 



consistently describes the full evolution of a rclativistic 
heavy-ion collision from the initial hard scatterings and 
string formation through the dynamical deconfinement 
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FIG. 14: (color online) The inverse slope parameter Te// of 
the dimuon yield from In+In at 158 A-GeV as a function of 
the dimuon invariant mass in PHSD compared to the data of 
the NA60 Collaboration [i,[z3- 



while simultaneously accounting for the electromagnetic 
radiation of the strongly coupled quark-gluon plasma 
(sQGP) via ofF-shell quark-antiquark annihilation, quark 
annihilation with gluon Bremsstrahlung and the gluon- 
Compton scattering mechanisms. 

In particular, the spectra in the intermediate mass 
range (1 GeV < M < 2.5 GeV) are found to be domi- 
nated by quark-antiquark annihilation in the nonpertur- 
bative QGP. Also, the observed softening of the trans- 
verse mass spectra at intermediate masses (1 GeV < 
M < 2.5 GeV) is approximately reproduced. 

Furthermore, for dileptons of low masses (M < 
0.6 GeV), a sizable contribution of partonic processes 
(in particular, the quark annihilation with the gluon 
bremsstrahlung) is found, thus possibly providing an- 
other window for probing the properties of the sQGP. 

Our present findings will have to be controlled by 
dilepton measurements at RHIC and LHC energies, since 
the PHSD approach is designed to operate also at these 
higher energies. Our results/predictions for higher ener- 
gies will be presented in near future. 



phase transition to the quark-gluon plasma as well as 
hadronization and to the subsequent interactions in the 
hadronic phase. 

In the present work, we have studied the dilepton 
production in In-I-In collisions at 158 A-GeV within 
the PHSD off-shell transport approach including a col- 
lisional broadening of vector mesons, microscopic sec- 
ondary multi-meson channels and the strongly interact- 
ing QGP radiation, which is described by the interac- 
tions of dynamical quasiparticles in line with the degrees 
of freedom propagated in the transport approach. 

A comparison to the data of the NA60 Collabora- 
tion shows that the dilepton yield is well described by 
including the coUisional broadening of vector mesons. 
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